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Crack  Opening  Displacement  of  Fatigue  Cracks  in  Supersolvus  INIOO 

B.  B.  Bartha^  R.  John^  M.  Caton^  S.  Jha^  W.  J.  Porter" 

Universal  Teehnology  Corporation,  Dayton,  OH  45433 

US  Air  Force  Research  Laboratory,  Wright-Patterson  AFB,  OH  45431 

University  of  Dayton  Research  Institute,  Dayton,  OH  45431 

Abstract 

Fatigue  crack  growth  tests  of  Supersolvus  INIOO  were  conducted  at  650C.  Multiple 
fatigue  cracks  were  identified  during  testing.  The  Crack  Opening  Displacement  (COD) 
of  each  fatigue  crack  was  investigated  using  Digital  Image  Correlation  (DIC).  The 
results  were  compared  to  Finite  Element  Analysis  (FEA)  results.  While  the  model  shows 
relatively  good  accuracy  for  longer  fatigue  cracks,  the  model  deviates  from  the 
experimental  results  for  the  short  fatigue  cracks  due  to  the  effect  of  local  microstructure 
and  load  history  effects. 

Introduction 

Fatigue  variability  plays  a  critical  role  in  the  determination  of  the  total  life  in  fracture 
critical  components  [refl].  Engine  Rotor  Life  Extension  (ERLE)  and  Materials  Damage 
Prognosis  programs  under  the  Air  Force  have  examined  various  life  limiting  factors  and 
the  potential  to  extend  the  lives  of  in  service  components.  Various  material  specific 
mechanisms  have  been  found  to  contribute  to  fatigue  variability.  In  titanium,  fatigue 
variability  was  shown  to  be  due  to  competing  failure  mechanisms  [ref2].  A  divergent 
failure  mechanism  has  also  been  reported  for  Rene’  88DT  superalloy  [ref  3].  The 
initiation  site  size,  location  of  the  initiation  site  and  the  crystallographic  orientation  near 
initiation  sites  were  shown  to  significantly  influence  the  total  fatigue  life  of  Rene’  88DT. 
(Discuss  Supersolvus  INIOO  fatigue  data  results  to  be  published). 
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While  previous  studies  have  investigated  the  fatigue  variability  of  niekel-base  superalloys 
foeusing  on  the  loeation  of  eraek  location,  the  fracture  mechanics  of  a  subsequent  fatigue 
crack  also  plays  an  important  role  in  determining  fatigue  variability.  The  current  study 
focuses  on  the  fatigue  crack  behavior  variability  in  supersolvus  IN  100  due  to  the  local 
microstructure  in  the  vicinity  of  a  fatigue  crack.  The  COD  behavior  of  each  fatigue  crack 
is  measured  with  optical  microscopy  combined  with  DIC  techniques.  The  COD  results 
are  used  to  compare  fatigue  cracks  that  propagated  at  a  different  rate  throughout  their 
fatigue  life,  while  having  the  same  nucleation  life. 

Include  DIC  technique  background  in  introduction  or  move  discuss  it  in  procedure? 

Experimental  Procedure 

A  fatigue  test  was  conducted  on  each  specimen  at  a  650C  temperature,  a  frequency  of 
0.33Hz,  a  stress  ratio  of  0.05  with  a  peak  stress  at  1  lOOMPa.  The  fatigue  experiment  was 
interrupted  at  cycle  numbers  1500,  2100,  2400,  2750  cycles  and  replicas  were  taken  of 
each  specimen  surface  at  room  temperature  to  identify  the  location  and  the  size  of  each 
fatigue  crack. 

The  specimen  surface  was  prepared  for  Electron  Backscatter  Diffraction  (EBSD)  and 
deformation  mapping  after  the  completion  of  the  fatigue  experiment.  Each  sample  was 
finish  polished,  and  fiducial  markers  were  placed  on  the  surface  of  each  specimen  (Figure 
xl).  The  5  micron  fiducial  marks  were  created  with  a  Vickers  indenter  at  a  50g  load 
setting.  EBSD  was  performed  in  the  area  of  each  crack  to  capture  the  grain  shape  and 
orientation  in  the  region  of  interest  at  a  1  micron  spot  size  resolution.  Each  sample  was 
etched  with  nitol?  solution  to  obtain  the  optimum  high  quality  optical  contrast  needed  for 
Digital  Image  Correlation. 

Deformation  mapping  was  performed  on  the  fatigue  cracks  of  interest  after  the 
completion  of  the  fatigue  test  (Figure  x2).  A  Nikon®  microscope  coupled  to  an 
Optronics®  digital  camera  was  used  to  capture  2048  by  2048  pixel  image  of  each  crack. 
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The  speeimen  was  loaded  and  unloaded  at  room  temperature  in  lOOMPa  stress 
inerements  to  1  lOOMPa.  Images  were  taken  of  each  crack  for  the  22  load  increments. 

Image  Analysis 

Digital  Image  Correlation  was  performed  on  each  image  with  ViC-2D  to  calculate  the 
displacement  field  surrounding  the  crack  for  each  load  increment.  The  image  at  the 
OMPa  unloaded  condition  was  used  as  a  reference  image  for  the  displacement  calculation. 
A  subset  size  of  99  pixels  was  used  for  each  calculation  with  a  5  pixes  step  size  and  cubic 
B-spline  interpolation.  The  resulting  displacement  data  was  combined  with  the  EBSD 
data  with  the  aid  of  the  fiducial  markers  to  align  the  optical  data  with  the  SEM  data. 

The  Crack  Opening  Displacement  (COD)  was  measured  from  the  displacement  data 
obtained  from  the  two  sides  of  the  crack  mouth.  The  displacements  calculated  are 
relative  displacements  from  the  zero  load  condition,  where  the  initial  position  of  the  crack 
mouth  is  unknown  from  the  experimental  results.  Therefore,  there  is  uncertainty  in  the 
total  COD  of  each  crack,  due  to  the  initial  unknown  COD  condition  that  cannot  be 
measured  experimentally. 

Results 

The  COD  measurements  were  compared  to  Linear  Elastic  Fracture  Mechanics  (LEFM) 
solution  from  Green  and  Sneddon  (Figure  x3)  [ref]  (Substitute  Franc  3D  results  for 
Greeen  Sneddon).  The  results  are  tabulated  in  Table  xl .  The  LEFM  solution  under¬ 
predicted  the  experimental  results  for  each  crack  (Figure  x4).  There  are  several  reasons 
for  the  deviation  of  the  LEFM  results  from  the  experimental  results.  The  LEFM  results 
do  not  account  for  the  three  dimensional  nature  of  the  problem  relating  to  specimen 
geometry,  and  the  local  microstructure.  Local  plasticity  effects  are  not  included  in  the 
calculations  near  the  crack  tip  and  from  the  first  shakedown  cycle  during  the  experiment. 
In  addition  the  COD  measurements  were  taken  at  room  temperature,  where  the  testing 
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was  conducted  at  elevated  temperature.  None  of  these  effects  are  accounted  for  in  the 
LEFM  calculations. 

The  experimental  results  were  also  compared  to  3D  FEA  results  using 
ABAQUS/ANSYS?  Coupled  with  Franc-3D.  The  results  show  that  the  FEA  predicted 
the  COD  values  more  accurately  for  the  larger  cracks  than  the  smaller  cracks.  The  COD 
was  significantly  over  predicted  by  the  3D  FEA  for  the  short  crack  cases  for  cracks  992E 
and  992E.  The  effect  of  local  microstructure  and  the  effect  of  closure  due  to  plasticity 
effects  are  considered  to  be  the  leading  causes  for  the  difference  in  the  COD  values  for 
the  short  cracks. 

Conclusion 

A  novel  optical  DIC  technique  was  used  to  measure  the  deformation  and  COD  of  fatigue 
cracks  in  supersolvus  IN  100.  The  results  were  compared  to  LEFM  and  3D  FEA  results 
to  show  the  significant  difference  in  the  COD  behavior  of  short  cracks  compared  to  long 
cracks.  The  difference  can  be  attributed  to  the  effect  of  local  grain  structure  and 
plasticity  in  the  vicinity  of  fatigue  cracks. 
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Figures  and  Tables 


Table  xl:  Comparison  of  measured  versus  theoretical  crack  size  and  the  number  of  load 
increments  taken  for  each  experiment. 


Crack 

Crack 
size  (lam) 

Exp 

(jim) 

GS 

(^m) 

Error 

(%) 

Inc 

(Mpa) 

Franc 

3D()xm 

Error 

(%) 

992A 

380 

1.48 

1.11 

25.0 

100 

1.48 

992C 

267 

1.03 

0.78 

23.9 

1100 

1.04 

-1 .5 

992D 

221 

0.96 

0.69 

28.1 

1100 

0.86 

10.4 

992E 

151 

0.73 

0.45 

39.0 

100 

0.59 

19.2 

979B 

366 

1.40 

1.07 

23.6 

100 

1.43 

-2.1 

979C 

341 

1.20 

1.00 

16.7 

1.33 

■EEl 

Figure  xl :  Optical  image  of  fatigue  crack  polished  and  etched  for  DIC  and  marked  for 
combinine  disnlacement  and  EBSD  results. 
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Figure  x2:  The  experimental  setup  eonsists  of  an  optieal  microseope  and  a  digital  camera 
coupled  and  mounted  to  a  four  axis  stage  focusing  on  a  region  of  interest  on  the  surface 
of  a  fatigue  specimen. 
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Q5-979  Crack  B 


05-979  Crack  C 
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Figure  x3:  Relative  COD  is  shown  at  1  lOOMPa  versus  0  MPa 
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05-992  Crack  A 


05-979  Crack  B 
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05-979  Crack  C 
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Figure  x4:  Peak  relative  crack  opening  displacement  for  loading  and  unloading  in 
lOOMPa  increments 


05-992  Crack  A 
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05-979  Crack  B 
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05-979  Crack  C 
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05-979  Crack  E 


Figure  x5:  Differential  COD  is  shown  relative  to  a  linear  fit  of  the  COD  slope  from 
hOOMPato  llOOMPa 
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05-992  Crack  A 
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05-979  Crack  B 
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05-979  Crack  C 
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05-979  Crack  E 
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Crack  Size  2a  (mm) 

Figure  x4:  The  COD  predietions  using  LEFM  are  consistently  lower  than  the  obtained 
experimental  values. 


19 


Chapter  2 


On  the  F/A-22  Full  Scale  Durability  Test  Results  Regarding  EIFS  Values 


The  caleulated  peak  stress  (99-155ksi)  in  the  analysis  of  the  cracks  in  TAR  45  are  near,  or 
over  the  yield  stress  (120ksi)  for  the  material.  The  shakedown  during  the  initial  plastic 
deformation,  and  strain  hardening  due  to  these  stresses  affects  not  only  the  EIFS,  but  also 
the  initiation  life  and  the  subsequent  crack  growth  properties  of  the  material.  The  test 
program  will  include  both  smooth  bar  and  notch  fatigue  tests  to  investigate  such  effects. 
10  uniaxial  smooth  bar  and  10  notch  fatigue  tests  (3  grain  sizes  ranging  from  0.015  up  to 
0.05  and  3  repeats  -  1  test  for  tuning  of  replication  procedures)  will  investigate  the 
experimental  EIFS,  possible  grain  boundary  and  size  effects,  as  well  as  initiation  and 
crack  growth  properties.  The  loading  range  and  notch  radius  will  be  chosen  to  produce 
stress  magnitudes  and  gradients  near  the  values  predicted  by  the  durability  testing. 

Acetate  replication  will  monitor  the  crack  initiation  and  size  at  various  points  during 
fatigue  life.  The  plastic  zone  size  will  be  determined  with  2-D  strain  mapping 
techniques. 

The  overview  of  Crack  2  in  TAR  45  mentions  that  the  determination  of  the  initiation  is 
obtained  from  the  local  stress  distribution  near  a  notch.  “Calculate  initiation  times  by 
multiplying  reference  stress  values  by  the  Boeing  shallow  gradient  (SSG)  factor  before 
determining  initiation  life.  The  stress  gradient  is  defined  as  the  delta  stress  in  the  first 
0.05  in  of  distance  from  the  peak  stress  divided  by  0.05.  ”  It  is  ambiguous  whether  the 
SSG  includes  previous  S/N  or  experimental  data.  It  is  also  unclear  how  the  spectrum 
loading  is  included  in  the  initiation  life  calculations.  Potentially,  a  comparison  between 
the  initiation  data  generated  from  the  above  fatigue  tests  and  the  initiation  life  analysis 
will  enhance  the  accuracy  of  the  initiation  life  prediction. 

The  peak  stress  near  the  location  of  the  cracks  is  predicted  by  a  coarse  FE  mesh  that 
brings  uncertainty  to  the  accuracy  the  peak  surface,  and  subsurface  stress  distribution. 
Unaccounted  shape  factor  effects  in  some  of  the  redesigned  components,  as  mentioned  in 
TAR  45,  further  decrease  the  accuracy  of  the  FE  model.  “The  lower  lug  fillet  was 
changed  to  remove  most  of  the  stress  concentration  that  caused  crack  2  at  A/C4008. 

There  is  no  detailed  solid  element  FEM  of  this  particular  configuration  and  so  the 
correlated  margins  and  lives  quoted  are  taken  from  BA3 12/Crack  1  to  be  conservative.  ” 
Higher  order  elements  may  not  be  sufficient  to  obtain  the  stress  distribution  accurately 
for  notches  that  have  high  peak  stress,  and  sharp  stress  gradients  with  respect  to  the 
overall  element  size.  A  convergence  study  of  the  peak  stress  and  the  stress  distribution  is 
necessary  to  determine  the  degree  of  mesh  refinement  needed  in  these  regions.  Sub¬ 
models  that  include  the  hotspots  are  necessary  for  verification.  The  accurate  subsurface 
stress  distributions  will  enhance  the  precision  of  both  the  initiation  and  the  crack  growth 
analysis. 
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Notch  Plasticity  effects  in  Supersolvus  INIOO 

B.  B.  Bartha^  R.  John^  M.  Caton^  S.  Jha^  W.  J.  Porter" 

Universal  Teehnology  Corporation,  Dayton,  OH  45433 

US  Air  Force  Research  Laboratory,  Wright-Patterson  AFB,  OH  45431 

University  of  Dayton  Research  Institute,  Dayton,  OH  45431 

Abstract 

A  tensile  test  of  Supersolvus  INIOO  was  conducted  at  room  temperature.  Multiple 
notches  were  placed  on  the  surface  of  the  specimen  using  a  Focused  Ion  Beam  (FIB). 
Deformation  mapping  was  conducted  on  each  notch  during  the  tensile  tests.  Digital 
Image  Correlation  (DIC)  was  performed  on  the  obtained  images  to  obtain  the  plasticity 
near  each  notch  area.  The  strain  results  were  compared  to  the  grain  orientation  data. 

Introduction 

Fatigue  variability  plays  a  critical  role  in  the  determination  of  the  total  life  in  fracture 
critical  components  [refl].  Engine  Rotor  Life  Extension  (ERLE)  and  Materials  Damage 
Prognosis  programs  under  the  Air  Force  have  examined  various  life  limiting  factors  and 
the  potential  to  extend  the  lives  of  in  service  components.  Various  material  specific 
mechanisms  have  been  found  to  contribute  to  fatigue  variability.  In  titanium,  fatigue 
variability  was  shown  to  be  due  to  competing  failure  mechanisms  [ref2].  A  divergent 
failure  mechanism  has  also  been  reported  for  Rene’  88DT  superalloy  [ref  3].  The 
initiation  site  size,  location  of  the  initiation  site  and  the  crystallographic  orientation  near 
initiation  sites  were  shown  to  significantly  influence  the  total  fatigue  life  of  Rene’  88DT. 
(Discuss  Supersolvus  INIOO  fatigue  data  results  to  be  published). 

While  previous  studies  have  investigated  the  fatigue  variability  of  nickel-base  superalloys 
focusing  on  the  location  of  crack  location,  the  fracture  mechanics  of  a  subsequent  fatigue 
crack  also  plays  an  important  role  in  determining  fatigue  variability.  The  current  study 
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focuses  on  the  crystal  plasticity  behavior  in  supersolvus  INI 00  due  to  the  local 
microstructure  in  the  vicinity  of  fatigue  notches.  The  crystal  plasticity  behavior  of  each 
FIB  notch  is  measured  with  optical  microscopy  combined  with  DIC  techniques.  The 
results  are  used  to  compare  the  FIB  notches  to  the  local  grain  structure  and  orientation. 

Experimental  Procedure 

A  tensile  test  was  conducted  on  a  specimen  at  a  room  temperature  to  a  peak  stress  at 
llOOMPa. 

Before  the  tensile  test,  he  specimen  surface  was  prepared  for  FIB  notch  placement,  and 
Electron  Backscatter  Diffraction  (EBSD).  Each  sample  was  finish  polished,  and  seven  40 
pm  FIB  notches  were  placed  150  mm  apart  on  the  surface  of  each  specimen.  EBSD  was 
performed  in  the  area  of  each  FIB  notch  to  capture  the  grain  shape  and  orientation  in  the 
region  of  interest  at  a  1  micron  spot  size  resolution.  The  sample  was  etched  with  nitol 
solution  to  obtain  the  optimum  high  quality  optical  contrast  needed  for  DIC 

Deformation  mapping  was  performed  on  the  FIB  notches  after  the  completion  of  the 
fatigue  test  (Figure  x2).  A  Nikon®  microscope  coupled  to  an  Optronics®  digital  camera 
was  used  to  capture  2048  by  2048  pixel  image  of  each  crack.  The  specimen  was  loaded 
and  unloaded  at  room  temperature  in  lOOMPa  stress  increments  to  1  lOOMPa.  Images 
were  taken  of  each  crack  for  the  22  load  increments. 

Image  Analysis 

Digital  Image  Correlation  was  performed  on  each  image  with  ViC-2D  to  calculate  the 
displacement  field  surrounding  the  crack  for  each  load  increment.  The  image  at  the 
OMPa  unloaded  condition  was  used  as  a  reference  image  for  the  displacement  calculation. 
A  subset  size  of  99  pixels  was  used  for  each  calculation  with  a  5  pixes  step  size  and  cubic 
B-spline  interpolation.  The  resulting  displacement  data  was  compared  to  the  EBSD  data. 

Results 
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Figures  1-14  show  the  local  strain  behavior  of  supersolvus  INIOO  around  seven  FIB 
notches.  While  the  nominal  strain  behavior  of  each  of  the  seven  areas  of  interest  were 
similar,  there  were  variations  in  the  local  strain  behavior  due  to  the  microstructure  near 
each  FIB  notch.  The  highest  local  strain  results  may  be  an  indicator  as  to  which  FIB 
notch  will  cause  a  crack  to  initiate  first.  Further  investigation  and  fatigue  testing  is 
necessary  in  this  case  to  obtain  the  fatigue  crack  growth  results,  and  to  compare  the  final 
fatigue  data  to  the  initial  plastic  behavior  of  the  seven  regions  for  the  first  cycle  tensile 
test. 
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Figure  1 :  Notch  1  strain  distribution  at  900MPa 
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Figure  2:  Notch  1  inverse  pole  figure  of  local  microstructure. 


Figure  3:  Notch  2  strain  distribution  at  900MPa 
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Figure  5:  Noteh  3  strain  distribution  at  900MPa 
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Figure  6:  Notch  3  inverse  pole  figure  of  local  microstructure. 
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Figure  7:  Notch  4  strain  distribution  at  900MPa 
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Figure  8:  Notch  4  inverse  pole  figure  of  local  microstructure. 
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Figure  9:  Notch  5  strain  distribution  at  900MPa 
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Figure  10:  Notch  5  inverse  pole  figure  of  local  microstructure. 
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Figure  1 1 :  Notch  6  strain  distribution  at  900MPa 
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Figure  12:  Notch  6  inverse  pole  figure  of  local  microstructure. 


Figure  13:  Notch  7  strain  distribution  at  900MPa 
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Figure  14:  Notch  8  inverse  pole  figure  of  local  microstructure. 
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In-Situ  Tensile  testing  of  Rene  88DT  at  room  and  elevated  temperatures 

B.  B.  Bartha^^,  S.  Fairchild'^,  T.  Murray‘S,  M.  Tschopp,  W.  J.  Porter‘s 
Universal  Teehnology  Corporation,  Dayton,  OH  45433 
US  Air  Force  Research  Laboratory,  Wright-Patterson  AFB,  OH  45431 
University  of  Dayton  Research  Institute,  Dayton,  OH  45431 

Abstract 

A  tensile  test  of  Rene  88  DT  was  conducted  at  70C  and  593C  inside  a  Scanning  Electron 
Microscope  (SEM).  Deformation  mapping  was  conducted  on  the  surface  of  the  specimen 
during  each  tensile  test.  Digital  Image  Correlation  (DIC)  was  performed  on  the  obtained 
images  to  obtain  the  plasticity  near  each  notch  area. 

Introduction 

Fatigue  variability  plays  a  critical  role  in  the  determination  of  the  total  life  in  fracture 
critical  components  [refl].  Engine  Rotor  Life  Extension  (ERLE)  and  Materials  Damage 
Prognosis  programs  under  the  Air  Force  have  examined  various  life  limiting  factors  and 
the  potential  to  extend  the  lives  of  in  service  components.  Various  material  specific 
mechanisms  have  been  found  to  contribute  to  fatigue  variability.  In  titanium,  fatigue 
variability  was  shown  to  be  due  to  competing  failure  mechanisms  [ref2].  A  divergent 
failure  mechanism  has  also  been  reported  for  Rene’  88DT  superalloy  [ref  3].  The 
initiation  site  size,  location  of  the  initiation  site  and  the  crystallographic  orientation  near 
initiation  sites  were  shown  to  significantly  influence  the  total  fatigue  life  of  Rene’  88DT. 
(Discuss  Supersolvus  IN  100  fatigue  data  results  to  be  published). 

While  previous  studies  have  investigated  the  fatigue  variability  of  nickel-base  superalloys 
focusing  on  the  location  of  crack  location,  the  fracture  mechanics  of  a  subsequent  fatigue 
crack  also  plays  an  important  role  in  determining  fatigue  variability.  The  current  study 
focuses  on  the  crystal  plasticity  behavior  in  Rene  88DT  due  to  the  local  microstructure  in 
at  room  and  elevated  temperatures. 
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Experimental  Procedure 

A  novel  in-situ  testing  teehnique  was  developed  to  obtain  the  loeal  deformation  behavior 
of  Rene  88DT  at  room  and  elevated  temperatures.  The  Experimental  setup  consisted  of  a 
Quanta  600  PEG  SEM  coupled  with  a  Fullam  lOOOlbs  high  temperature  in-situ  loading 
stage  (Figure  1).  A  specimen  was  prepared  with  a  gage  length  was  measured  to  be  0.1 1  in 
wide  and  0.032in  thick  with  a  0.7in  gage  length  (Figure  2).  Through  thin  film  ablation 
TTFA  was  used  to  coat  the  specimen  surface  with  platinum  nano-particles,  in  order  to 
obtain  images  with  enough  contrast  to  conduct  image  correlation  on  the  obtained 
secondary  electron  images  (Figure  3-4). 

The  specimen  surfaces  were  cleaned  to  a  finish  1  pm  mechanical  polish  or  electro- 
polished  before  the  platinum  deposition.  It  was  vital  in  these  experiments  to  keep  the 
surface  clean  of  everything  for  optimum  imaging  and  deposition.  The  surface 
microstructure  can  be  seen  in  the  electron  backstatter  micrograph  in  Figure  5.  TTFA 
consisted  a  a  lOnm  platinum  target  deposited  onto  a  fused  silica  plate  that  was  transparent 
to  the  laser  wavelength  of  248nm  with  an  energy  density  of  0.5  J/cm^2.  The  chamber 
was  filled  with  argon  at  a  pressure  of  5  torr. 

Each  tensile  specimen  was  mounted  into  the  Fullam  stage  grip  fixture  and  preloaded  to 
1  Olbs  to  allow  the  specimen  to  seat  and  to  mitigate  specimen  rotation.  The  images  were 
collected  at  16  bit  depth  at  a  pixel  resolution  of  4096  by  3773  pixels.  The  stage  controller 
was  turned  off  during  the  image  acquisition  of  each  load  step  to  minimize  distortion 
effects  crated  by  the  motor  operating  the  stage.  Each  image  was  focused  by  raising  and 
lowering  the  stage  fixture  in  order  to  keep  the  working  distance  the  same.  If  focus  is 
performed  with  the  beam  only,  rotational  artifacts  have  been  detected  with  DIG. 


Image  Analysis 
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Digital  Image  Correlation  was  performed  on  each  image  with  ViC-2D  to  calculate  the 
displacement  field  for  each  load  increment.  The  image  at  the  OMPa  unloaded  condition 
was  used  as  a  reference  image  for  the  displacement  calculation.  A  subset  size  of  99 
pixels  was  used  for  each  calculation  with  a  5  pixels  step  size  and  cubic  B-spline 
interpolation. 

Results 

The  displacements  and  strains  were  calculated  with  each  load  step  and  can  be  seen  in 
Figure  6.  The  strain  concentrations  correlate  well  to  the  observed  slip  bands  that  can  be 
seen  in  Figure  7.  The  slip  bands  in  the  region  of  interest  were  activated  at  different  load 
levels,  where  the  slip  planes  need  to  be  obtained  from  the  Electron  Backscatter  data  seen 
in  Figure  8.  A  Schmid  factor  analysis  is  needed  to  correlate  the  Schmid  factor  values  to 
the  local  strain  concentrations  seen  from  DIG. 

Current  plasticity  models  compare  closely  with  the  obtained  results,  and  can  be  used  to 
improve  their  accuracy.  These  results  are  vital  in  improving  the  current  crystal  plasticity 
models  in  this  area. 


Figure  1 :  High  temperature  in-situ  tensile  stage 
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Figure  3:  Seeondary  electron  micrograph  of  platinum  deposit  on  Rene  88  DT  surface. 


Figure  4:  Secondary  electron  micrograph  of  the  platinum  coating  used  for  DIC 
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S  Microscope  Accelerating  Voltage  Working  Distance  Detector 

Quanta  PEG  15kV  12.6  mm  SSD  — 50  pm 


Figure  5:  Backseatter  electron  micrograph  showing  the  local  grain  structure 
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Figure  6:  Local  strain  distribution  at  986MPa  calculated  from 
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Figure  7:  Seeondary  electron  micrograph  showing  the  slip  bands  formed 
At  986MPa 
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MICROSTRUCTURAL  CONDITIONS  CONTRIBUTING  TO  FATIGUE  VARIABILITY 
IN  P/M  NICKEL-BASE  SUPERALLOYS 
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Abstract 

Variability  in  the  fatigue  behavior  of  two  eommon  niekel-base 
superalloys  is  diseussed,  with  emphasis  given  to  understanding 
the  behavior  of  short  fatigue  eraeks.  Rene88DT  and  INI 00  are 
the  materials  of  interest  and  serve  important  roles  in  many  turbine 
engine  systems.  Multiple  speeimens  of  eaeh  material  were  tested 
under  low  eyele  fatigue  (LCF)  eonditions  at  elevated  temperature. 
Two  INI 00  speeimens  that  exhibited  signifieant  differenee  in 
eyeles  to  failure  and  three  Rene88DT  speeimens  having  dissimilar 
short  fatigue  eraek  growth  rates  were  interrogated  to  determine 
the  reasons  for  these  variations.  This  paper  will  diseuss  the 
mierostrueture  adjaeent  to  the  initiation  sites  in  these  speeimens 
and  its  role  in  affeeting  the  observed  disparate  eraek  growth 
behavior.  Additionally,  an  image  eorrelation  teehnique,  used  to 
evaluate  the  development  of  strain  eoneentrations  in  a  Rene88DT 
tensile  speeimen  at  room  temperature,  is  deseribed.  The  utility  of 
image  eorrelation  for  identifying  mierostruetural  ‘hot  spots’  (i.e. 
initiation  sites)  in  eoneert  with  LCF  testing  is  eonsidered. 

Introduction 

Life  predietion  models  that  integrate  the  inherent  fatigue 
variability  of  materials  used  in  gas  turbine  engines  are  of  interest 
to  the  USAF.  Past  researeh  indieates  that  mueh  of  the  variability 
seen  in  fatigue  behavior  is  related  to  material  response  in  the 
presenee  of  small  (<  100pm)  fatigue  eraeks  [1,2].  While  extensive 
databases  of  long  eraek  growth  data  exist  for  niekel-base 
superalloys,  similar  databases  for  small  eraek  growth  data  are 
laeking,  as  is  our  general  understanding  of  their  behavior 
partieularly  under  LCF  eonditions.  Faetors  assoeiated  with 
differenees  in  small  eraek  behavior  inelude  loeal  mierostruetural 
features  ineluding  neighboring  grain  (mis)orientation,  grain 
substrueture  (y’  preeipitates),  twins,  and  anomalies  sueh  as  the 
presenee  of  inelusions  and  pores  [3-7]. 

Rene88DT  and  INI 00  are  important  alloys  in  the  field  of  Ni-base 
powder  metal  superalloys  in  turbine  engine  applieations.  They  are 
the  foeus  of  this  paper  due  to  their  important  role  in  US  Air  Foree 
(USAF)  propulsion  systems.  A  eurrent  foeus  of  USAF  materials 
researeh  pertains  to  the  development  of  life  predietion 


methodologies  based  on  determining  remnant  eomponent 
eapability  via  prognostie  approaehes.  The  ultimate  goal  of  this 
work  is  to  arrive  at  a  better  understanding  of  those  faetors  eentral 
to  influeneing  fatigue  variability  to  aid  future  USAF  life- 
management  efforts. 

Experimental 

Material 

The  alloys  investigated  in  this  study  were  Rene88DT  and  IN  100. 
Both  are  powder  metal  proeessed  materials  that  were  solution 
treated  above  the  gamma  prime  (y’)  solvus  temperature  (super- 
solvus).  The  sourees  of  material  for  eaeh  alloy  were  large 
paneake  forgings  produeed  using  standard  industrial  praetiees. 
The  nominal  eomposition  and  mierostrueture  of  eaeh  alloy  are 
shown  in  Table  I  and  Figure  1,  respeetively.  The  average  grain 
size  of  the  Rene88DT  of  approximately  30  pm  (Fig.  la);  the 
mierostrueture  is  highlighted  by  a  large  number  of  annealing 
twins  and  very  fine  y’  phase  material  that  forms  upon  eooling  with 
diameters  on  the  order  of  200  nm  (Fig.  lb).  Aging  or  tertiary  g’ 
phase,  although  not  apparent  in  Fig  lb,  is  also  present  and 
reported  to  be  approximately  lOOA  in  size  [8]. 

The  grain  size  of  the  INI 00  is  approximately  25  pm  (Fig  le).  A 
bi-modal  distribution  of  y’  partieles  is  formed  upon  eooling  from 
annealing,  with  the  eooling  y’  having  diameters  of  1pm  (Fig.  Id) 
and  the  aging  y’  having  diameters  of  ~  lOOnm.  Twinning  is 
apparent  in  the  INI 00,  but  not  as  prevalent  as  in  the  Rene88DT. 

The  tensile  properties  eorresponding  to  the  temperatures  at  whieh 
the  fatigue  testing  was  eondueted  for  eaeh  alloy  follow.  At  593C, 
the  Rene88DT  has  a  0.2%  offset  yield  of  985  MPa  and  an  elastie 
modulus  of  176  GPa  [3].  At  650C,  the  INI 00  measured  a  0.2% 
yield  of  888  MPa,  ultimate  tensile  stress  of  1225  MPa,  an  elastie 
modulus  of  175  GPa,  and  a  strain  at  failure  of  41%. 

Testing  Proeedure 

Cylindrieal  dog-bone  speeimens  having  a  nominal  gage  diameter 


Table  1.  Nominal  eomposition  of  Rene88DT  and  INI 00  (Weight  Pereent) 


Co 

Cr 

Mo 

W 

A1 

Ti 

V 

N 

B 

c 

Zr 

Ni 

R88DT 

13 

16 

4 

4 

2.1 

3.7 

0.  7 

0.01 

0.03 

0.03 

bal 

INI  00 

18.5 

12.4 

3.2 

5.5 

4.7 

0.8 

0.07 

bal 
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of  6.35  mm  were  used  in  this  study.  The  speeimens  were 
extraeted  sueh  that  the  tensile  axes  were  parallel  to  the 
eireumferenee  of  eaeh  forging.  All  of  the  speeimens  were 
maehined  to  a  low  stress  grind  (LSG)  finish.  The  INI 00 
speeimens  were  tested  in  the  LSG  surfaee  eondition.  To  further 
minimize  the  effeets  of  surfaee  residual  stresses  due  to  maehining, 
the  gage  seetions  of  the  ReneSSDT  speeimens  were 
eleetropolished  using  a  proeedure  deseribed  in  [3]. 

The  ReneSSDT  speeimens  were  tested  under  load  eontrol  at 
1250MPa,  593°C,  a  frequeney  of  0.33  Hz,  and  a  stress  ratio  of 
0.05.  All  INI 00  testing  was  also  done  under  load  eontrol  at 
llOOMPa,  650°C,  a  frequeney  of  0.33  Hz  and  a  stress  ratio  of 
0.05.  The  stress  levels  ehosen  for  eaeh  group  of  speeimens  are 
representative  of  low  eyele  fatigue  (LCF)  eonditions. 

Replieate  testing  of  ReneSSDT  speeimens  at  1250  MPa  using  the 
eonditions  previously  deseribed  in  eoneert  with  the 
eleetropolished  gage  seetions  resulted  in  a  predominanee  of 
surfaee-nueleated  eraeking.  The  surfaee-nueleated  eraeking 


served  as  a  eonvenient  platform  for  investigating  variable  short 
eraek  growth  rates  using  aeetate  replieation  teehniques.  Surfaee 
replieas  were  taken  periodieally  (at  least  every  thousand  eyeles, 
often  more  frequently)  throughout  the  tests  to  monitor  eraek 
advanee.  Replieation  was  done  at  room  temperature  and,  as  sueh, 
required  the  fumaee  to  be  shut  off  and  the  sample  eooled  to 
ambient. 

An  array  of  analytieal  tools  was  used  to  examine  the  role  of 
mierostrueture  on  the  initiation  and  propagation  of  small  fatigue 
eraeks.  These  tools  ineluded:  1)  a  seanning  eleetron  mieroseope 
(SEM)  for  fraetography  and  orientation  imaging  mieroseopy 
(OIM)  of  material  at  and  adjaeent  to  the  initiation  sites,  2)  a 
foeused  ion  beam  (FIB)-equipped  SEM  to  reveal  the 
mierostrueture  in  the  areas  of  interest,  3)  a  software  paekage  that 
enabled  quantitative  fraetographie  measurements  of  the  initiation 
sites  and  surrounding  material  (MeX  by  Alieona  Imaging  GmbH) 
[9]  and  4)  a  eharge-eoupled  deviee  (CCD)  eamera  attaehed  to  an 
optieal  mieroseope  for  digital  image  eorrelation  (DIC). 


Figure  1.  Baekseattered  SEM  images  of  a,  b)  ReneSSDT  and  e,d)  IN  100  mierostruetures 
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Rene88DT  Fatigue 

Surface  replicas  were  used  to  monitor  the  advance  of  short 
cracks  in  three  Rene88DT  specimens.  The  initiation  site  for 
each  of  these  was  found  to  be  non-metallic  inclusions  (NMI) 
located  at  the  specimen  surface.  The  lengths  of  each  NMI 
exposed  at  the  specimen  surface  are  shown  in  Table  2.  The 
NMI  were  similar  in  size,  shape  and  composition.  Fracture 
surfaces  for  each  specimen  are  shown  in  Figure  2.  The  distinct 
thumbnail  appearance  and  banding  around  each  initiation  site 
seems  to  be  an  artifact  of  the  replication  process  which  required 
periodic  cool  down  to  room  temperature.  The  bands  can 
highlight  the  amount  of  crack  advance  during  a  given  cycle 
block.  However,  as  shown  in  Figure  3,  for  the  specimen 
identified  as  ‘slow’,  there  were  blocks  that  showed  signs  of 
crack  arrest  or  negligible  growth  (e.g.  from  2000  to  6000 
cycles). 


Table  2:  Information  for  cracks  in  Rene88DT  selected 
for  characterization. 


Crack  ID 

Initiation 

Type 

Surface  length 
of  initiating 

feature  (|dm) 

Cycles  to 
120(im 

Slow 

NMI 

35 

-7700 

Medium 

NMI 

45 

-4500 

Fast 

NMI 

40 

-2100 

In  an  effort  to  correlate  fracture  surface  roughness  at  crack 
lengths  (2a)  less  than  200  pm  to  crack  growth  rate,  stereopairs 
of  each  fracture  surface  were  collected  for  input  to  the  MeX 
software.  A  three-dimensional  representation  (digital  elevation 
map)  of  each  surface  was  created  and  area-based  surface 
roughness  measurements  were  made.  A  summary  of  the  surface 
roughness  measurements  are  shown  in  Table  3.  Using  a  line- 
intercept  technique,  the  projected  area  is  the  planar  measurement 
of  the  area  of  interest  (i.e.  length  x  width).  An  example  of  a  grid 
placed  over  the  initiation  area  of  the  ‘slow’  specimen  is  shown 
in  Figure  4.  The  true  area  is  a  summation  measurement  of  the 
changes  in  topography  along  the  lines  of  the  grid  and  root  mean 
square  (RMS)  height  is  a  value  calculated  by  taking  the  square 
of  the  difference  in  height  from  a  mean  height.  The  qualitative 
statement,  “all  else  being  equal  (i.e.  material,  testing  conditions, 
etc.),  the  more  tortuous  a  fracture  surface,  the  longer  the  life”  is 
supported  by  the  quantitative  findings  reported  in  Table  3.  The 
trends  shown  in  Table  3,  the  higher  the  true/  projected  ratio  and 
the  higher  the  RMS  value,  the  slower  the  crack  will  advance 
hold  true. 


Table  3.  Summary  of  surface  roughness  measurements  near  the 
_ initiation  sites  for  Rene88DT  specimens. _ 


projected 
area  (pm^) 

true  area 
(pm^) 

Ratio  of 
true/ 

projected 

RMS 

height 

(pm) 

Slow 

19200 

22557 

1.175 

3.109 

Medium 

18186 

20697 

1.138 

2.271 

Fast 

19248 

21433 

1.114 

2.093 

As  mentioned  earlier,  variable  growth  rates  exhibited  by  small 
fatigue  cracks  are  often  attributed  to  neighboring  grain 


20pni 

b) 


c) 


Figure  2.  Fracture  surfaces  and  initiation  sites  from  Rene88DT 
specimens,  a)  slow,  b)  medium,  c)  fast.  Note  arrows  on  b) 
indicating  grain  boundary. 
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misorientations.  A  good  example  of  this  is  shown  in  Figure 
2b.  The  arrows  on  Figure  2b  indieate  the  point  at  whieh  two 
river  lines  terminate,  indieating  the  presenee  of  a  grain  boundary 
between  two  grains  with  appreeiable  misorientation.  Craek 
advanee  is  signifieantly  slowed  at  this  point  in  turn  hindering 
movement  of  the  entire  eraek  front.  The  eraek  growth  rate  slows 
in  this  region  as  shown  in  Figure  3  for  the  ‘medium’  speeimen 
from  approximately  90  to  100  pm  (eraek  length,  2a,  at  the 
surfaee).  This  findng  is  in  agreement  with  loeation  of  the 
arrows  in  Figure  3b  whieh  are  found  to  be  approximately  45pm 
(eraek  length,  a,  for  subsurfaee)  from  the  initiation  site. 
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Figure  3.  Surfaee  eraek  size  vs  N  for  ReneSSDT. 


Figure  4.  Grid  used  for  line-intereept  measurement  of  fraeture 
surfaee  roughness  in  ReneSSDT  ‘slow’  speeimen. 


fNlOO  Fatigue 

Twenty  one  (21)  IN  100  speeimens  were  tested  in  this  portion  of 
the  study  with  fatal  eraeks  initiating  at  either  sites  of  porosity 
(70%)  or  non-metallie  inelusions  (30%).  Most  of  the  fatal 
eraeks  initiated  at  the  surfaee  of  the  speeimen  (60%).  due  to  the 
large  number  of  pores  as  eompared  to  NMI,  the  vast  majority  of 
the  surfaee  initiated  eraeks  were  assoeiated  with  porosity.  From 
this  group,  two  speeimens  that  initiated  at  pores  on  the  speeimen 
surfaee  were  seleeted  for  eharaeterization.  While  the  pore 
morphology  and  loeation  within  eaeh  speeimen  were  very 
similar,  a  roughly  4X  differenee  in  eyeles  to  failure  (5105  vs 
23,233)  was  exhibited.  Information  about  the  pores  from  eaeh 
speeimen  is  shown  in  Table  4.  One  point  of  interest  is  that  the 
speeimen  initiating  at  the  larger  pore  had  longer  life. 


Table  4:  Size  and  morphology  of  pores  measured  in  INI  00 
_ speeimens  seleeted  for  eharaeterization.  _ 


Speeimen 

ID 

Cyeles  to 
failure 

Pore 

Diameter 

(4m) 

Pore 

Area 

(4m 

Aspeet 

ratio 

long  life 

23,233 

32.5 

S31 

1.1 

short  life 

5,105 

20.3 

326 

1.0 

The  next  step  was  to  use  the  FIB  to  prepare  the  material  beneath 
and  surrounding  the  initiation  sites  for  OIM  imaging.  Figure  5b 
and  e  show  post-FIB,  baekseattered  eleetron  images  of  eaeh  of 
the  IN  100  speeimens  revealing  the  supersolvus  mierostrueture. 
Figure  5e  and  f  show  the  eleetron  baekseattered  diffraetion 
(EBSD)  images  for  eaeh  speeimen. 

As  expeeted,  given  the  P/M  proeessing  route  used  for  the  IN  100, 
the  OIM  information  from  eaeh  speeimen  reveals  a 
mierostrueture  with  no  strong  overall  texture.  The  most 
noteworthy  differenee  in  the  OIM  data  between  the  two 
speeimens  is  the  presenee  of  more  twin  boundaries  in  the  long 
life  speeimen.  Efforts  are  eontinuing  to  determine  whether  this 
differenee  in  the  amount  of  twinning  in  the  short  eraek  region  is 
a  faetor  eontributing  to  variable  fatigue  behavior.  Other  faetors 
under  investigation  inelude  quantitative  analysis  using  the  OIM 
data  of  the  types  of  grain  boundaries  present  in  the  areas  of 
interest.  A  preponderanee  of  high  angle  (large  misorientation) 
versus  low  angle  (small  misorientation)  grain  boundaries  has 
been  shown  to  substantially  hinder  eraek  growth  rates  and  may 
play  a  signifieant  role  in  fatigue  variability. 

ReneSSDT  Tensile/  Image  Correlation 

In  addition  to  the  fatigue  variability  study  with  mierostrueture, 
tensile  testing  was  performed  while  imaging  the  deformation 
behavior  of  the  loeal  mierostrueture.  Eleetron  Baekseatter 
Diffraetion  (EBSD)  was  first  performed  in  an  area  of  interest  on 
the  surfaee  of  a  ReneSSDT  speeimen  in  order  to  obtain  the  loeal 
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100|jm 


d) 


100|jm 


c)  f) 

Figure  5.  SEM  images  of  the  IN  100  fraeture  surfaees  before  FIB  milling,  (a,  d),  after  FIB  milling,  (b,  e),  and  OIM  images  from  the 
FIB-milled  surfaees,  (e,  f).  Figs.  a,b,&e)  eorrespond  to  the  long-life  speeimen;  Figs.  d,e,&  f  to  the  short-life  speeimen. 


orientation  (Figure  6).  The  orientation  values  were  used  to  obtain 
the  grain  boundary  misorientation  (Figure  7).  During  the  tensile 
test,  images  were  taken  of  the  area  of  interest  with  a  CCD  eamera 
and  an  optieal  mieroseope.  The  strain  distribution  ealeulated 
using  Vie-2D  digital  image  eorrelation  software  (Figure  8)  [10, 
11] 

.  The  EBSD  data  was  eombined  with  the  optieal  data  to  eompare 
the  effeet  of  orientation  on  the  plastie  strains.  The  highest  strain 
eoneentrations  are  loeated  near  the  boundary  of  the  large  grain 
shown  in  the  image.  While  several  twins  within  the  large  grain 
show  elevated  strain  levels,  most  of  the  high  strain  areas  oeeur 
near  high  misorientation  angle  grain  boundaries. 

The  grain  boundary  misorientation  direetly  affeets  the  miero- 
meehanieal  behavior  in  the  above  mierostrueture.  Strain 
eoneentrations  near  twiee  the  mean  strain  or  higher  are  observed 
near  high  angle  misorientation  grain  boundaries.  Other  faetors 


sueh  as  neighboring  grain  shapes  and  sizes  and  orientation  with 
respeet  to  the  loading  direetion  also  eontribute  to  the  observed 
elevated  strain  levels. 

Due  to  experimental  limitations,  the  data  presented  is  on  the 
surfaee  only  two  dimensional,  even  though  the  problem  is  three 
dimensional.  Therefore,  the  subsurfaee  grain  properties  and  the 
out  of  plane  displaeements  are  unknown.  While  the  ability  to 
measure  these  three  dimensional  properties  would  inerease  the 
aeeuraey  and  the  eonfidenee  in  the  experimental  results,  the  trends 
due  to  grain  boundary  misorientation  and  twinning  would  stay  the 
same.  The  twinning  also  produeed  some  strain  eoneentrations  of 
signifieanee  when  eompared  to  grain  boundary  misorientations. 
Further  investigation  at  higher  magnifieations  is  neeessary  to 
eompare  and  eontrast  the  loeal  strain  behavior  of  grain  boundaries 
and  twin  boundaries  and  their  importanee  in  eraek  initiation  and 
propagation  behavior  relating  to  fatigue  variability. 
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Conclusions 


Figure  6:  Inverse  pole  figure  of  area  imaged  during  deformation 
mapping. 


—  «•  «• 
50gm 

Figure  7:  Grain  boundary  misorientation  map  near  large  eenter 


Figure  8:  Total  strain  field  of  area  of  interest  with  a  mean  total 
strain  at  10%  at  1400MPa  tensile  stress. 


Rene88DT  and  INI 00  LCF  speeimens  that  exhibited  dissimilar 
small  fatigue  eraek  growth  rates  and  variable  eyeles  to  failure, 
respeetively,  were  examined  to  determine  the  mierostruetural 
eonditions  responsible  for  eaeh  behavior.  A  digital  image 
eorrelation  teehnique  was  employed  on  a  Rene88DT  tensile 
speeimen  to  determine  how  and  where  plastie  strains  develop 
during  loading. 
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Abstract 

Electron  backscatter  diffraction  (EBSD)  data  was  obtained  for  the  primary  and  secondary 
alpha  phases  of  Ti-6Al-2Sn-4Zr-6Mo  (Ti-6246).  The  results  established  that  the 
observed  microtexture  in  duplex  alloys  is  a  direct  result  of  the  prior  beta  grain 
orientations,  and  variant  selection.  In  addition,  for  a  homogeneous  duplex 
microstructure,  all  of  the  secondary  alpha  phase  and  the  majority  of  the  primary  alpha 
grains  retain  crystallographic  coherency,  according  to  the  Burgers  Relationship,  with  the 
surrounding  beta-phase  matrix  which  comprise  prior  beta  grains.  This  investigation 
develops  a  unique  technique  to  separate  the  phases  in  a  duplex  a/p  titanium 
microstructure.  The  image  quality  map  is  used  to  separate  the  two  phases  based  on  the 
gray  level,  size  and  aspect  ratio  for  each  phase. 
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Introduction 

The  study  of  mierotexture  in  duplex  titanium  alloys  with  a  microstructure  that 
consists  of  equiaxed  (primary)  and  platelet  (secondary)  alpha  phase  that  is  surrounded  by 
a  thin  layer  of  beta  phase  at  ambient  temperatures  is  complicated  by  several  factors. 
First,  in  order  to  examine  if  there  exists  a  correlation  in  the  crystallographic  orientation  of 
the  alpha  phase  constituents  (primary  and  secondary  alpha)  with  that  of  the  beta  phase, 
high  resolution  EBSD  patterns  are  required  to  collect  a  statistically  representative  number 
of  data  from  each  phase.  Second,  in  some  alloys,  such  as  Ti-6A1-4V,  where  the  volume 
fraction  of  retained  beta  phase  at  ambient  temperatures  is  very  low,  it  is  even  more 
difficult  to  measure  the  crystallographic  orientation  of  the  beta  phase. 

From  a  mechanical  perspective,  mierotexture  is  an  artifact  of  thermomechanical 
processing  that  can  have  a  deleterious  effect  on  the  fatigue  properties  of  titanium  alloys. 
Recent  work  examining  the  role  of  crystallographic  orientation  of  the  primary  alpha 
phase  on  fatigue  behavior  has  shown  that  primary  alpha  oriented  suitably  for  basal  <a> 
slip  serves  as  the  fatigue  crack  initiation  sites  [1].  Invariably,  the  microstructural  volume 
surrounding  the  fatigue  crack  initiation  site  was  found  to  be  strongly  textured  for  basal 
and  prism  slip.  Hence  from  a  fatigue  life  prediction  perspective,  tools  that  aid  in  the 
understanding  of  how  microtextures  develop  and  the  ability  to  predict  crystallographic 
orientation  are  essential  for  the  development  of  life  prediction  models  that  incorporate 
crystallographic  orientation  into  a  modeling  framework. 
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The  objective  of  the  present  work  is  to  use  image  quality  maps  to  separate  the 
primary  alpha  phase  and  the  secondary  alpha  phase  and  study  the  texture  of  the  two  phase 
microstructure. 

Results  and  Discussion 

Figures  1-5  show  the  equiaxed  and  lath  alpha  phase  obtained  from  EBSD  data  Each  set 
of  a  set  of  4X4  EBSD  data  in  Figure  6  was  processed  to  produce  only  a  set  of  grains  with 
only  equiaxed  or  lath  alpha,  and  the  prior  beta  grains  were  predicted  using  the  previous 
Monte  Carlo  Method.  Figures  7-9  show  the  predicted  prior  beta  grains  from  the  original 
microstructure,  as  well  as  the  equiaxed  and  lath  microstructures  only.  The  lath 
microstructure  clearly  shows  a  large  prior  beta  grain  in  red.  While  the  equiaxed 
microstructure  only  shows  a  fraction  of  the  coherency  of  the  previous  data  set.  The 
results  show  that  The  equiaxed  alpha  does  retain  some  coherency  to  the  prior  beta  phase 
of  the  microstructure. 


Figure  1 :  Inverse  pole  figure  showing  the  grain  orientation  of  a  duplex  microstructure. 
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Figure  2:  Obtain  Euler  angle  for  eaeh  primary  alpha  grain 


Figure  3:  Equiaxed  alpha  grains  grown  from  original  mierostrueture 
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Figure  4:  Fill  in  Euler  angles  of  secondary  alpha  phase  from  the  original  EBSD  data. 


Figure  5:  Lath  alpha  grains  grown  from  original  microstructure. 
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Figure  6:  Composite  inverse  pole  figure  map  of  a  4X4  grid  of  EBSD  Data 


Figure  7:  Transformed  beta  of  the  entire  data  set 
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Figure  8:  Transformed  beta  of  the  lath  alpha  grains  only 


Figure  9:  Transformed  beta  of  the  equiaxed  alpha  grains  only. 
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Abstract 

Electron  backscatter  diffraction  (EBSD)  data  was  obtained  for  the  primary  and  secondary 
alpha  phases  of  Ti-6Al-2Sn-4Zr-6Mo  (Ti-6246).  The  results  established  that  the 
observed  microtexture  in  duplex  alloys  is  a  direct  result  of  the  prior  beta  grain 
orientations,  and  variant  selection.  In  addition,  for  a  homogeneous  duplex 
microstructure,  all  of  the  secondary  alpha  phase  and  the  majority  of  the  primary  alpha 
grains  retain  crystallographic  coherency,  according  to  the  Burgers  Relationship,  with  the 
surrounding  beta-phase  matrix  which  comprise  prior  beta  grains.  This  investigation 
develops  a  unique  technique  to  separate  the  phases  in  a  duplex  a/p  titanium 
microstructure.  The  image  quality  map  is  used  to  separate  the  two  phases  based  on  the 
gray  level,  size  and  aspect  ratio  for  each  phase. 
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Introduction 

The  study  of  mierotexture  in  duplex  titanium  alloys  with  a  microstructure  that 
consists  of  equiaxed  (primary)  and  platelet  (secondary)  alpha  phase  that  is  surrounded  by 
a  thin  layer  of  beta  phase  at  ambient  temperatures  is  complicated  by  several  factors. 
First,  in  order  to  examine  if  there  exists  a  correlation  in  the  crystallographic  orientation  of 
the  alpha  phase  constituents  (primary  and  secondary  alpha)  with  that  of  the  beta  phase, 
high  resolution  EBSD  patterns  are  required  to  collect  a  statistically  representative  number 
of  data  from  each  phase.  Second,  in  some  alloys,  such  as  Ti-6A1-4V,  where  the  volume 
fraction  of  retained  beta  phase  at  ambient  temperatures  is  very  low,  it  is  even  more 
difficult  to  measure  the  crystallographic  orientation  of  the  beta  phase. 

From  a  mechanical  perspective,  mierotexture  is  an  artifact  of  thermomechanical 
processing  that  can  have  a  deleterious  effect  on  the  fatigue  properties  of  titanium  alloys. 
Recent  work  examining  the  role  of  crystallographic  orientation  of  the  primary  alpha 
phase  on  fatigue  behavior  has  shown  that  primary  alpha  oriented  suitably  for  basal  <a> 
slip  serves  as  the  fatigue  crack  initiation  sites  [1].  Invariably,  the  microstructural  volume 
surrounding  the  fatigue  crack  initiation  site  was  found  to  be  strongly  textured  for  basal 
and  prism  slip.  Hence  from  a  fatigue  life  prediction  perspective,  tools  that  aid  in  the 
understanding  of  how  microtextures  develop  and  the  ability  to  predict  crystallographic 
orientation  are  essential  for  the  development  of  life  prediction  models  that  incorporate 
crystallographic  orientation  into  a  modeling  framework. 
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The  objective  of  the  present  work  is  to  use  image  quality  maps  to  separate  the 
primary  alpha  phase  and  the  secondary  alpha  phase  and  study  the  texture  of  the  two  phase 
microstructure. 

Results  and  Discussion 

Figures  1-14  show  the  step  by  step  numerical  procedure  that  was  implemented  on  each 
IQ  map  image  in  order  to  separate  the  primary  and  secondary  alpha  phase  and  the  prior 
beta  phase  using  the  gray  scale  information.  The  contrast  was  initially  stretched  in  order 
to  find  the  gray  scale  level  that  separates  the  primary  alpha  phase  from  the  beta  phase. 

The  subsequent  areas  were  binned  and  marked  with  a  grain  ID,  and  the  size  and  aspect 
ratio  of  each  grain  was  calculated.  The  results  were  used  to  further  filter  out  smaller 
grains  as  well  as  secondary  alpha  phase  grains  with  a  high  aspect  ratio.  The  resulting 
image  contained  only  the  primary  alpha  phase.  The  Euler  angle  for  each  grain  was 
obtained  from  the  original  EBSD  scan. 

The  secondary  alpha  phase  was  found  in  a  similar  manner.  Initially  the  gray  level  from 
the  primary  alpha  phase  was  separated  from  the  original  IQ  map  image.  The  image 
contrast  was  again  changed  to  accurately  subtract  the  beta  phase.  The  smaller  grains 
were  again  subtracted  out,  and  grains  with  low  aspect  ratios  were  also  ignored  to  find  the 
secondary  beta  phase.  The  Euler  angles  for  the  resulting  grains  were  obtained  from  the 
original  EBSD  data. 

Figures  15-17  show  the  pole  figures  of  the  total  image,  the  primary  alpha  phase,  and  the 
secondary  alpha  phase  that  was  obtained  from  the  cleaned  data  set.  The  difference  in  the 
texture  in  the  different  phases  near  crack  initiation  areas  is  important  to  characterize  in 
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order  to  understand  failure  mechanisms  and  to  develop  accurate  life  prediction 
methodologies  for  duplex  microstructures. 


Figure  1 :  IQ  Map  of  duplex  Microstructure. 


Figure  2:  Contrast  Stretch  o  IQ  map 
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Figure  3:  Threshold  grayscale  to  obtain  alpha  phase 


Figure  4:  Subtract  out  small  grains 
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Figure  5:  Subtract  out  grains  with  high  aspect  ratios 


Figure  6:  Resulting  primary  alpha  microstructure  after  image  processing 
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Figure  7:  Obtain  Euler  angle  for  eaeh  primary  alpha  grain 


Figure  8:  Subtract  gray  scale  values  from  original  IQ  map 
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Figure  9:  Stretch  out  the  contrast  of  each  image 


Figure  10:  Threshold  image  to  subtract  beta  phase. 
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Figure  1 1 :  Subtract  out  small  grains 


Figure  12:  Subtract  out  grains  with  low  aspect  ratios 
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Figure  13:  Resulting  seeondary  alpha  phase  grains. 


Figure  14:  Fill  in  Euler  angles  of  seeondary  alpha  phase  from  the  original  EBSD  data. 
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Figure  15:  Pole  figure  of  the  entire  data  set. 
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Figure  16:  Pole  figure  of  the  Primary  alpha  phase  only. 
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Figure  17:  Pole  figure  of  the  seeondary  alpha  phase  only. 
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1.  Introduction 

Sandwich  panels  are  used  extensively  in  aircraft  structures  due  to  the  very  low  weight  to 
stiffness  ratios  produced  using  various  materials  and  geometries.  Metallic,  ceramic  and 
composite  materials  are  each  popular  in  specific  loading  and  thermal  conditions  where 
they  are  the  most  beneficial  structurally,  and  cost  effectively.  The  application  of 
sandwich  panels  to  Thermal  Protection  Systems  (TPS)  for  hypersonic  flight  is  currently 
an  area  of  interest  where  the  use  of  metallic  tiles  may  be  advantageous  due  to  the  high 
damage  tolerance  of  the  material  over  ceramics  and  composites  that  are  in  use  for  TPS 
systems.  The  performance  of  the  metallic  tiles  must  be  determined  and  compared  to 
existing  TPS,  and  the  cost  effectiveness  of  metallic  panels  needs  to  be  investigated  to  find 
if  metallic  TPS  is  a  viable  option  over  current  ceramic  and  composite  materials. 

Various  modeling  and  characterization  techniques  on  metallic  TPS  panels  are  needed  to 
determine  their  performance  and  mechanical  behavior.  The  results,  in  turn,  can  be  used 
to  determine  the  optimum  weight  and  shape  of  these  panels.  Extensive  research  is 
available  covering  structural  modeling  using  methods  from  the  theory  of  elasticity  as  well 
as  Finite  Element  Methods  (FEM).  The  modeling  efforts  include  structural  models  based 
on  sandwich  property  and  geometry,  loading  configurations  based  on  mechanical  and 
thermal  loads  and  attachment  boundary  condition  issues.  Dynamic  issues  such  as 
vibrations,  fatigue,  and  impact  are  also  areas  of  interest  for  modelers  and  designers.  The 
results  from  these  models  aid  designers  choose  the  correct  material,  panel  shape,  for  a 
particular  attachment  and  mechanical  and  thermal  load  state. 

In  order  to  develop  structural  models  for  weight  optimization  of  honeycomb  panels,  it  is 
necessary  to  understand  currently  available  modeling  techniques  and  conditions  that  drive 
models  as  they  relate  to  metallic  TPS  panels.  The  relevant  modeling  and  characterization 
techniques  can  be  chosen,  modified  and  further  developed  for  weight  optimization  of 
metallic  TPS  panels. 

2.  Background  on  Modeling  Techniques 

Modeling  of  sandwich  panels  has  been  developed  using  theory  of  elasticity  with  beam 
theory.  The  equations  from  beam  theory  give  accurate  description  of  the  displacement 
and  stress  fields  of  sandwich  and  honeycomb  structures  for  various  boundary  conditions 
and  load  cases.  The  results  have  been  used  to  analyze  bending,  buckling  and  wrinkling 
failure  criteria  [1].  In  addition,  anisotropy,  dynamic  effects,  as  well  as  thermal  and 
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moisture  effects  have  been  developed  to  solve  more  complex  issues  for  composite 
materials  and  added  boundary  conditions  [2],  Layered  and  non-linear  analysis  have  also 
been  applied  to  anisotropic  layered  composites  [3].  Most  models  for  specific  applications 
apply  specific  cases  of  these  theories,  incorporating  specific  conditions  that  are  important 
for  each  application  that  the  basic  theories  may  not  account  for  accurately.  These 
advances  in  the  field  have  been  reviewed  periodically  over  the  years  [4],  The  overview 
of  the  different  sections  of  the  field  is  presented  here  to  emphasize  the  issues  that  arise  in 
honeycomb  structures,  and  the  possible  role  that  they  play  in  weight  optimization 
problems. 

Static  Loads 

Honeycomb  structures  have  superior  strength  under  bending  and  compressive  load 
configurations  because  they  exhibit  their  largest  strength  in  out  of  plane  loading. 

Buckling  and  bending  are  common  modes  of  failure  in  these  loading  configurations.  The 
structural  behavior  of  these  modes  have  been  modeled  using  plates  and  shells  in  FEM  as 
well  as  approximated  analytical  solutions  [5,  6].  Several  failure  modes  such  as  yield, 
shear  or  wrinkling  from  3-Point  bending  of  the  honeycomb  structure  can  also  occur  for 
various  effective  density  and  core  thickness  configurations  [7].  In  addition,  plasticity 
solutions  have  also  been  developed  in  conjunction  with  elastic  models  to  predict  the 
macroscopic  deformation  behavior  of  sandwich  panels  incorporating  kinematic  hardening 
and  softening  [8].  Many  of  the  analytical  models  have  been  validated  experimentally 
using  three  and  four  point  bend  fixtures  to  a  high  degree  of  accuracy  [9].  The  relative 
density  of  the  honeycomb  panels  also  directly  correlates  to  the  buckling  stiffness  [10,  11]. 

Ultimate  stress  conditions  are  also  applied  to  existing  elastic  models  to  investigate  and 
address  shear,  buckling,  debonding  and  fracture  failure  modes  [12].  When  3-D  triaxial 
loads  are  applied  to  honeycomb  structures,  various  modes  of  failure  dominate  depending 
on  the  shape  and  the  orthotropy  of  the  sandwich  structure  [13].  Asymmetry  of  the 
cellular  structures  has  a  significant  effect  on  the  strength  of  structures  that  is  tied  to 
apparent  density  the  composite  [14].  More  complex  curved  sandwich  panels  have  also 
been  used  for  structural  leading  edge  applications.  Computational  models  are  used  for 
panels  that  these  exhibit  unique  3-D  properties  [15]. 

Dynamic  Loads 

Honeycomb  panels  may  undergo  cyclic  loading  that  can  cause  failure  in  the  panels  or 
fastener  attachments.  In  some  of  these  cases,  the  stress  field  is  found  with  beam  solutions 
while  approximating  the  fatigue  properties  using  Paris  law  for  crack  propagation, 

Basquin  law  for  high  cycle  fatigue  and  Coffm-Manson  law  for  low  cycle  fatigue  [16]. 
Other  properties  of  sandwich  panels  relating  to  fatigue  strength  of  fasteners  attachments 
have  also  been  investigated  experimentally  [17].  The  inherent  mechanical  behavior  of  the 
material  used  for  honeycomb  panels  in  these  cases  needs  to  be  determined  before  the  use 
of  the  above  modeling  techniques. 
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Modeling  the  load  history  of  honeyeomb  panels  subjeeted  to  impaet  loading  is  needed  to 
determine  the  damage  toleranee  of  the  panels.  Considerable  work  has  been  conducted  on 
the  vibrational  behavior  of  sandwich  panels  from  impact  to  dynamic  response.  Both 
experimental  and  analytical  studies  have  been  conducted  investigating  the  free  vibration 
of  honeycomb  panels  [18].  There  have  also  been  studies  on  the  strength  of  various 
honeycomb  panels  subjected  to  impulse  loading  and  designs  of  optimum  sandwich  cores 
for  shock  resistance  [19].  The  vibrational  response  of  these  materials  can  also  be 
modeled  using  Timoshenko  beam  theory,  as  well  as  relating  the  response  to  the  material 
elastic  constants  [20]. 

Honeycomb  structures  are  also  used  in  thermal  applications  such  as  heat  exchangers.  In 
these  cases,  it  is  important  to  understand  the  thermomechanical  characteristics  of  the 
material  in  addition  to  the  strength  of  the  structures  because  there  is  usually  an  inherent 
tradeoff  between  the  performance  of  each  behavior  [21].  Large  thermal  gradients  can 
also  cause  large  deformations  and  thermal  strains  on  panels  that  lead  to  failure.  In  these 
cases  the  value  of  the  Coefficient  of  Thermal  Expansion  (CTE),  and  structural  heat 
capacity  must  be  chosen  carefully  for  each  material  in  order  to  mitigate  out  of  tolerance 
deformations,  and  to  avoid  debonding  issues  that  are  associated  with  composite  materials. 

Materials 

Composite  materials,  as  well  as  layered  and  smart  structures  such  as  piezoelectric 
honeycombs  are  heterogeneous  and  difficult  to  model,  and  need  modifications  to 
analytical  methods  to  solve  for  the  displacements  and  stress  fields  [22].  Homogenization 
methods  have  been  used  in  these  cases  for  unit  cell  composites  [23,  24].  These  models 
convert  the  heterogeneous  properties  of  the  composite  honeycomb  structure  to  a 
homogeneous  effective  elasticity  tensor  and  plate  properties.  Laminate  plate  theory 
assumptions  are  used  in  applications  where  multiple  layers  and  materials  are  used  to 
build  honeycomb  panels  [25,  26].  Strain  energy  and  weight  function  methods  can  also  be 
used  to  model  composite  structure  behavior  along  with  the  homogenization  methods  for 
single  unit  cells  [27]. 

Metallic  sandwich  structures  are  advantageous  because  of  their  strength  to  weight  ratios 
are  combined  with  their  added  damage  tolerance  over  composite  and  ceramic  structures 
[28].  Improved  strength  in  metallic  honeycomb  structures  for  bending  and  buckling 
resistance  can  also  be  gained  with  fdler  foam  cores  instead  of  only  thickening  column 
walls  [29].  The  combined  strength  and  thermal  characteristics  are  advantageous  for 
specific  applications  where  the  structure  must  have  multi-purpose  use  [30]. 

TPS  Applications 

Metallic  TPS  have  shown  to  have  a  promising  combination  of  thermal  and  structural 
performance  characteristics  that  are  applicable  in  space  vehicle  applications  [31,  32]. 
Weight  optimization  of  metallic  TPS  is  an  important  factor  when  comparing  existing 
ceramics  and  composites  [33].  In  addition,  the  thermal  performance  of  the  metallic  TPS 
must  be  considered  making  sure  the  material  does  not  operate  above  critical 
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temperatures,  and  has  proper  CTE,  heat  eapacity,  and  heat  flux  charaeteristies.  Thermal 
bending  and  heat  shielding  eaused  by  the  high  surfaee  temperatures  and  subsequent  steep 
subsurfaee  temperature  gradient  are  issues  that  must  be  eonsidered  in  TPS  designs.  In 
addition,  attachment  issues  relating  to  the  performance  of  the  panel  attachment  region, 
and  the  cost  effectiveness  of  replacing  TPS  panels  with  particular  attachments  need  to  be 
resolved. 

A  combination  of  the  above  modeling  techniques  can  also  be  used  to  develop  weight 
optimization  models  to  obtain  the  correct  material  and  sandwich  configurations  for 
several  loading  cases  [1,  34].  Many  parameters  can  be  varied  such  as  honeycomb 
geometry,  while  others  such  as  the  material  properties  can  be  fixed  to  optimize  the 
problem  [35,  36,  37].  Studies  have  shown  that  analytical  solutions  for  optimum  design 
can  be  obtained  if  the  constraints  are  up  to  three  parameters  for  certain  cases,  and 
numerical  analysis  is  needed  when  considering  additional  variables  [38]. 

3.  Weight  Optimization  Application 

The  development  of  a  comprehensive  weight  optimization  model  for  honeycomb 
sandwich  TPS  must  include  the  main  factors  that  drive  the  performance  of  the  panels.  It 
is  vital  to  the  accuracy  of  the  model  to  know  the  panel  geometry,  the  type  of  attachments, 
as  well  as  the  input  structural  and  thermal  loads  on  each  panel.  These  geometrical  and 
load  variables  define  the  boundary  conditions  for  both  elasticity  and  finite  element  based 
models.  The  knowledge  of  mechanical  properties  of  the  chosen  materials  is  needed  to 
predict  the  stress  state  of  the  panels  for  each  boundary  condition.  The  performance 
characteristics  are  needed  to  predict  the  failure  locations  for  a  particular  stress  state.  The 
combined  variables  of  these  configurations  and  material  properties  determine  the  panel 
material,  geometry,  and  weight  for  this  application. 

Initially,  elasticity  based  beam  theory  for  honeycomb  panels  will  be  used  to  optimize  the 
weight  of  the  honeycomb  panel  similar  to  previous  work  performed  for  TPS  [31].  The 
constraints  will  include  peak  stress  values  for  bending  and  buckling  failure  modes.  The 
optimization  model  will  also  need  to  account  for  attachment  boundary  conditions  by 
using  pinned,  fixed  or  free  boundary  conditions  that  change  the  peak  stress  level  for  a 
given  load.  Structural  and  thermal  loading  input  boundary  conditions  can  then  be 
included  in  the  model  by  including  point,  and  distributed  loads,  as  well  as  thermal  strain 
and  thermal  gradient  effects.  Variation  in  material  properties  will  need  to  be  included 
using  FEA  for  more  complex  models.  Commercial  software  package  such  as  ISight  will 
be  used  to  optimize  the  weight  for  a  particular  set  of  the  described  input  variables  using 
numerical  routines  such  as  the  downhill  simplex  method. 

Many  of  the  variables  used  in  the  described  weight  optimization  scheme  must  be 
constrained  to  obtain  a  feasible  solution.  The  results  of  the  research  efforts  from 
AFRL/ML,  AFRL/VA  and  industry  will  aid  in  choosing  if  variables  such  as  materials  and 
temperatures  can  be  constrained,  and  determining  input  bounds  such  as  web  and  panel 
thickness.  The  resulting  weight  optimization  model  will,  in  turn,  provide  a  tool  to 
determine  the  structural  feasibility  of  using  metallic  TPS  that  is  competitive  with 
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materials  used  in  current  applications.  The  model  will  act  as  a  design  tool  to  determine 
optimum  attachment  configurations,  and  material  usage  for  TPS  applications. 

Weight  Optimization  Trade  Study  With  Hypersizer 
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Figure  1 :  Unit  Weight  of  a  30  in  IN617  panel  with  various  cell  sizes  under  increasing 
normal  pressures. 
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Figure  2:  Unit  weight  of  a  30  in  IN617  panel  with  various  cell  sizes  under  increasing 
normal  pressures  for  a  0. 1  in  normal  constraint. 
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